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Abstract. The processing of lunar regolith for the production of oxygen is a key component of 
the In-Situ Resource Utilization plans currently being developed by NASA. In the carbothermal 
process, a portion of the surface of the regolith in a container is heated by exposure to a heat source 
so that a small zone of molten regolith is established. A continuous flow of methane is maintained 
over the molten regolith zone. In this paper, we discuss the development of a chemical conversion 
model of the carbothermal process to predict the rate of production of carbon monoxide. Our model 
is based on a mechanism where methane pyrolyzes when it comes in contact with the surface of the 
hot molten regolith to form solid carbon and hydrogen gas. Carbon is deposited on the surface of 
the melt, and hydrogen is released into the gas stream above the melt surface. We assume that the 
deposited carbon mixes in the molten regolith and reacts with metal oxides in a reduction reaction 
by which gaseous carbon monoxide is liberated. Carbon monoxide bubbles through the melt and 
is released into the gas stream. It is further processed downstream to ultimately produce oxygen. 

Keywords: Carbothermal processing, lunar regolith, pyrolysis, kinetics 
PACS: 47.70.Fw, 82.40.-g 
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NOMENCLATURE 


a 

A 

r 

v --ra 

k 

P 

P* 

r p 

t, t* 

t C 


ti 


t 


m 


to 


v 

V 0 

Xco 
a, a* 

P, P * 
Po 


7 

5 

A 

A' 


stoichiometric coefficient in chemical reaction in Eq (1) 
melt surface area 

concentration of methane (mole/m 3 ) 
reaction rate constant (m 3 / (mole • s)) 

variable proportional to the square root of time, p = (2kp 0 t*) 1 ^ 2 

value of p when methane flow is initiated, p^ = (2k/3 0 t i ) 1 ^ 2 

flux of carbon from pyrolysis (mole/ (m 2 • s)) 

scaled and physical (s) time, respectively 

scaled time at which CO production rate is a maximum 

time when methane flow is initiated (s) 

scaled time beyond which methane flow rate is turned off 

reference time for a growing melt (s) 

melt volume (m 3 ) 

reference melt volume for a growing melt (m 3 ) 
number of moles of carbon monoxide (mole) 

scaled and physical (mole/m 3 ) concentration of carbon in the melt, respectively 

scaled and physical (mole/m 3 ) concentration of metal oxide in the melt, respectively 

initial concentration of metal oxide in the melt (mole/m 3 ) 

rate constant for the pyrolysis of methane (m/s) 

scaled excess moles of carbon in the melt, 5 = A t m — 1 

dimensionless parameter, A = 

dimensionless parameter, A' = (2kPoto) l ^ 2 \ 


1 Introduction 


NASA’s plans for lunar exploration involves the production of oxygen for life support and use 
as a propellant. Oxygen is available in the lunar regolith as oxides of various materials. The 
extraction of gaseous oxygen from the oxides by reduction using hydrogen (Taylor and Carrier, 
1992), molten salt electrolysis (Tripuraneni-Kilby et al., 2006) and carbothermal processing of 
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regolith (ORBITEC, 2006) are currently being pursued by NASA (Sanders et al., 2008), and is a 
prime candidate for In-Situ Resource Utilization. 

Carbothermal processing refers to the use of carbon in some form, together with high processing 
temperature, to reduce the oxides. Examples of carbothermal processes on Earth are the production 
of silicon from silica using wood, charcoal or coal (a standard process), and the production of 
titanium-rich materials by refining low-grade ilmenites using graphite (Wang and Yuan, 2006). 

The carbothermal processing of lunar regolith using methane as the source for carbon was first 
performed by Rosenberg and co-workers at Aerojet in the 1960’s (Rosenberg et al., 1991). It 
has been revisited by researchers at Orbital Technologies Corporation (ORBITEC) since the early 
1990’s. An excellent reference and resource for the application of the carbothermal process for the 
conversion of lunar regolith to ultimately produce oxygen is a report by ORBITEC (ORBITEC, 
2006; see also Gustafson, White and Fidler, 2009). ORBITEC used lunar regolith simulants in 
their processing experiments. 

Figure 1 shows a schematic sketch of the carbothermal reactor. Lunar regolith is contained in a 
suitable container. A portion of the surface of the regolith is heated by exposure to a heat source 
such as a laser beam or a concentrated solar heat flux, so that a small zone of molten regolith is 
established within the container. The unmelted regolith surrounding the molten zone thermally 
insulates it from the container. A heat shield is used above the molten zone to minimize radiative 
heat loss. A continuous flow of methane is maintained over the molten regolith zone. Since 
the molten regolith is very hot (temperature will exceed 1650°C), methane will pyrolyze when it 
contacts the molten zone. Carbon is deposited on the surface of the melt, and hydrogen is released 
into the gas phase above the melt surface. The deposited carbon reacts with metal oxides (chiefly 
Si0 2 ) in the molten regolith in a reduction reaction in which carbon monoxide is liberated. CO 
bubbles through the melt and is released in the gas phase above the melt. It is carried downstream 
for further chemical processing to ultimately yield oxygen, which is liquified and stored. Our 
objective is to model the production of carbon monoxide in the carbothermal reactor, and predict 
the rate at which CO is liberated as a function of the various physicochemical parameters of the 
system. We also pursue simple strategies to control the buildup of carbon in the melt. 
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FIGURE 1: Sketch of the Carbothermal Reduction of Lunar Regolith. 

2 Model formulation and assumptions 

We assume that carbon deposition occurs on the hot molten surface via a surface reaction. This 
surface reaction is regarded as first order with respect to a gaseous carbon-containing species and 
produces solid carbon and gaseous hydrogen. This carbon-containing species need not necessarily 
be methane; it can be some other carboneous species which may be produced by the gas-phase 
pyrolysis of methane. In reality, depending on the operating pressure, flow conditions and tem- 
perature profiles above the hot molten surface, there would be a distribution of carbon-containing 
species (e.g., methane, ethylene, acetylene, methyl radical, etc.,) and each of them would react 
on the surface to deposit carbon at their own rates. Our model, however, “lumps” all of these 
possibilities into a “single” species and describes the deposition process by a first-order surface re- 
action of this “predominant” species. For sake of simplicity, we call this “lumped” species methane 
throughout the paper. Thus, the flux of carbon deposited on the surface of the melt is assumed to 
be proportional to the concentration of methane in the gas phase adjacent to the melt surface. 

We assume that the rate of transport of gaseous species from the bulk flow to the melt surface 
is not limiting the surface reaction rate. This is expected to be valid in the practical application, 
where the goal would be to maximize the carbon deposition rate so that the CO production rate 
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is maximized. This can be accomplished by an ample and continuous supply of methane in the 
bulk flow so that its rate of diffusion through the boundary layer above the melt surface is suffi- 
ciently fast to replenish its consumption at the surface. Under these conditions, the concentration 
of methane at the surface can be assumed to be constant and approximated by its value at the flow 
inlet. The surface reaction rate constant needs to be extracted by direct comparison of model pre- 
dictions with experimental data. The analysis can easily be extended to situations where gas-phase 
transport effects are significant. 

The deposited carbon is assumed to be ‘well-mixed’ in the melt such that there are no spatial 
gradients of the concentration of carbon or the metal oxide within the melt. This would be the 
case when the zone is sufficiently small such that diffusion of carbon within the melt is rapid rela- 
tive to the time of its reaction with the metal oxide. Mixing is more likely caused by temperature 
gradients in the melt that generate a flow due to natural convection and thermocapillary forces. 
Besides, the reaction in the melt liberates gaseous carbon monoxide in it, and the gravity-induced 
bubble motion is likely to render the melt to be well-mixed. The melt viscosity is expected to play 
a significant role in determining the relative rates of carbon advection and diffusion. Heiken et cil. 
(1991) mention that a synthetic liquid with the composition of a lunar basalt sampled by Apollo 1 1 
has a viscosity in the range 4.5 poise (1495°C) to 10 poise (1395°C) and that terrestrial lavas are 
ten times more viscous. The viscosity of melts typically depends on the processing temperature 
and composition, in particular the proportion of MgO, CaO, FeO and Fe 3 0 4 , etc., in comparison 
with those of Si0 2 and A1 2 0 3 (Themelis, 1995). 

The carbon in the melt is assumed to undergo a single-step, first-order chemical reaction with the 
metal oxide present in the molten regolith. This reaction rate constant also needs to be determined 
by comparison of the model predictions with experimental data. We assume that the deposited 
carbon does not form a crust on the surface of the melt, shielding the melt surface from the heat 
flux and the methane flow. While crusting has been observed under some conditions (ORBITEC 
2006), rapid mixing due to bubble evolution and controlled rate of deposition of carbon greatly 
diminishes the likelihood of crust formation. 
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We assume that the molten zone is established and maintained by heating the regolith by suit- 
able means, in a batch-processing mode. The heat transfer aspects of the problem are not analyzed 
here. The processes of chemical conversion and heat transfer within the molten regolith are as- 
sumed to be uncoupled. Let V and A denote the volume and the surface area of the molten zone 
of regolith. A is assumed to be constant and independent of time during the chemical conversion 
process. We will analyze cases where the melt volume V is (i) a constant independent of time, and 
(ii) increases as the square root of time. A constant melt volume implies that the molten regolith 
fully occupies a fixed volume of its container (the regolith shown in Figure 1 is fully molten and 
the container is in contact with the melt). It might also be applied for a case where the heat transfer 
process is such that V and A are slowly varying in time compared to the time scale for chemical 
conversion. When there is contact between the melt and the solid regolith particles, we assume 
that there is a well-defined interface demarcating the melt and the solid regolith particles. Thus 
no trickling of the melt into the regolith bed is considered. This assumption appears to be true in 
the experiments of Gustafson et al. (2009) where they report their post-experiment observation of 
only loose binding of the regolith particles to a solidified molten mass. 

The analysis is organized as follows. We first consider V to be constant, and determine the evolu- 
tion of the concentration of carbon and metal oxide in the melt. This is followed by an analysis of 
a possible strategy to minimize the amount of leftover carbon in the melt. Finally, a growing melt 
where V increases as the square root of time, for constant A, is also analyzed. 

Let a* and (5* denote the molar concentration (moles/m 3 ) of carbon and metal oxide, respec- 
tively, within the melt. At time t* = 0, the initial concentrations are a* = 0 and /3* = To- The 
chemical reaction between carbon and the metal oxide in the regolith is 

a MO x + C ^CO + a M. (1) 

The stoichiometric coefficient of the oxide term in the above reaction is a = l/x. We assume 
that the chemical reaction is first order with respect to carbon and the metal oxide, with a rate 
constant k (m 3 / (moles • s)). The flux of carbon due to pyrolysis of methane is denoted by r p , 
and is assumed to be linearly proportional to the molar concentration of methane C m in the gas 
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phase adjacent to the melt surface, with the constant of proportionality denoted by 7 (m/s). The 
equations for the evolution of a* and (3* in the melt can be obtained from mass conservation to be: 


da* 

dt* 


A^Cm 

V 


-ka* (3*, 


d/3* 

dt* 


—a k a* (3* . 


The instantaneous rate of production of carbon monoxide is 

dX co _ _Vdf3* 

dt* a dt* 


( 2 ) 

(3) 

(4) 


3 Results and Discussion 


3.1 Constant melt volume 


We assume that the melt area and volume are constant. The methane concentration C m in the gas 
phase is also assumed to be constant. Dimensionless variables are defined as follows. 


a* (3* t* 

(3 0 /a ^ /V f l/(k(3o) 


(5) 


The scaled equations and boundary conditions are 


da 

dt ^ 


(6) 


d P O 

M = - Ul3 


(7) 


a(0) = 0, 0 ( 0 ) = 1 (8) 

where A = ' s a dimensionless parameter. A is a measure of the relative importance of the 

rate of production of carbon by the surface pyrolysis of methane to its rate of consumption by 
reaction with the metal oxide. The former rate is constant for typical operating conditions of fixed 
temperature, pressure, and sufficient flow rate to maintain a constant methane concentration. The 
latter rate varies depending on the accumulation or depletion of carbon and the availability of the 
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metal oxide in the melt. So, the magnitude of A indicates the relative dominance of the two rates in 
a competitive process; i.e., small lambda means relatively small carbon deposition rates, and vice 
versa. 


Equation (4) for the instantaneous rate of production of carbon monoxide can be expressed as 

dX co Vkffi 


dt* 


a p. 


(9) 


Thus the instantaneous rate of production of carbon monoxide is proportional to the product of a 
and p. The solution can be determined to be 


P = 


(y, — f3 — l + Xt 


2a/X Exp (t — |t 2 ) 


2^ + Exp (i) Erf ( 4= ) + Erf 


tA-1 


( 10 ) 


( 11 ) 


V2\J 1 V V2A, 

The following expression is an asymptotic expansion for P that is especially useful for small values 
of A. 

1 - Xt 


P = 


(A - 3A 2 + . . .)(! - At)Exp(-t + ft 2 ) + (l - + (i^yr + • • •) 


( 12 ) 


(1-Ai) 2 ^ (1— Ai) 4 

Plots of the scaled concentrations of carbon (a) and metal oxide (P), and the scaled rate of pro- 
duction of carbon monoxide (which is proportional to the product a P, see Equation 9) from the 
solution given above are shown in Figure 2 for various values of A. It is worthy of note that, in the 
simplified model that we have developed, the only parameter that determines the behavior of the 
reactor is A. This parameter captures the effects of melt size relative to its surface area, reaction 
rate constant, pyrolysis rate constant, and the concentrations of methane and metal oxide on the 
rate at which chemical conversion takes place in the reactor. 


Figure 2 shows that the metal oxide concentration in the melt decreases monotonically with time, 
and the amount of carbon in the melt increases monotonically with time. The rate of production 
of carbon monoxide is initially zero and increases with time as the reaction between the deposited 
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(a) A = 0.01 




(b) A = 0.1 (c) A = 1 



(d) A = 10 

FIGURE 2: Plot of the scaled concentrations of carbon (a) and metal oxide ((3) in the melt and 
the scaled rate of production of carbon monoxide (a/3) versus scaled time (t) for various values of 
A. 
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(a) (b) 

FIGURE 3: (a) the scaled time t c at which the carbon monoxide production rate is a maximum and 
(b) the corresponding scaled maximum rate of production of CO (a(t c )/3(t c )). both as a function 
of A. 

carbon and the metal oxide progresses, and finally decreases to zero as the reaction is completed. 
The exact shape of the curves in Figure 2, especially that for the the rate of production of CO, is 
quite sensitive to the value of A. 

As the metal oxide is depleted, the reaction finally slows down and stops. Figure 2 shows that 
the rate of production of CO attains a maximum before it eventually decreases to zero. From 
Equations (10) and (11), it can be shown that the maximum is attained at a time t — t c which 
satisfies the following equation. 

O r -1 \ - 1 1/2 

/5(C) = i (l-Af c )+ [-(l-AC) 2 + -j (13) 

t c represents the time for which there is an optimal balance between the accumulation of carbon in 
the melt due to deposition at the melt surface and its consumption by reaction with the metal oxide. 
The dependence of t c and the scaled maximum CO production rate a(t c )/3(t c ) on A are shown in 
Figure 3. t c monotonically decreases with increasing A, while the scaled maximum CO production 
rate is a monotonically increasing function of A. The results for t c would be useful in design and 
operation of the reactor to ensure that requirements on the amount and rate of production of CO 
are met, and identify times beyond which there is diminishing returns in operating the reactor. 
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Since the amount of carbon in the melt increases with time, there is unreacted carbon that is left 
behind in the melt after the reaction is completed. A useful strategy to minimize the amount of 
unreacted carbon is to turn the supply of methane off at an appropriate time before the reaction is 
completed, such that the residual amount of carbon present in the melt at that time is sufficient to 
react with the remaining amount of the metal oxide. This scenario is addressed below. 

3.2 Control of carbon in the melt 

It can be seen from Figure 2 that towards the end of the conversion process (for large t), a is 
quite large. This represents the large amount of carbon that has been deposited on the surface of 
the melt. On the other hand, the concentration of metal oxide in the melt is quite low, showing 
that much of the metal oxide has been converted. Therefore, much of the carbon that has been 
deposited remains unreacted in the melt. Once the chemical conversion process is terminated after 
a sufficient amount of the metal oxide in the regolith is converted, the unreacted carbon remains 
in the spent regolith. It is wasted, unless the spent regolith is specifically processed to recover the 
carbon. Since carbon (in the form of methane) is envisioned to be transported from Earth to the 
moon, it is desirable to minimize the carbon loss by minimizing the unreacted carbon in the melt. 

One strategy to minimize the carbon loss is by turning the flow of methane off at some point in 
time before the end of the conversion process. The idea is that the amount of carbon present in the 
melt at the time the methane flow is turned off is sufficient to react with the residual metal oxide 
in the melt. Since there is no fresh carbon that is deposited to the melt after the flow of methane 
is turned off, the steep increase in the amount of carbon in the melt depicted in Figure 2 will be 
avoided. 

By turning the methane flow off, we mean that the concentration of methane in the inlet gas stream 
is zero. Since there are other sub-systems present downstream of the carbothermal reactor to 
process the carbon monoxide produced, there must still be gases flowing through the entire system. 
We anticipate that even with the concentration of methane in the inlet stream to the carbothermal 
reactor being zero, the gas flow rate through the entire system is maintained by employing other 
‘sweep’ gases (such as hydrogen, carbon monoxide, inert gas, etc.). 
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We assume that the methane flow is turned off at t — t 0 . For t < t 0 , the solution given in Equations 
(10) and (11) applies. The chemical conversion for t > t 0 can be modeled as 


da 

dt 


—a (3 


(14) 


d(3 

dt 


—a [3 


(15) 


a(t o) = a to , p(t 0 ) = (3 to 


(16) 


where 


Ao — 


Exp (t 0 - |fo) 


(17) 


2v /A + Exp (i) [Erf j + Erf 

Note that a to = (3 to — 1 + A t 0 . Let a to — (3 to = Xt 0 — l — S. 6 represents the moles of carbon in 
the melt in excess of stoichiometry at the time that the methane flow is turned off. The solution to 
Equations (14) - (16) is 

(3 t0 5 e S{t0 ~^ 


P = 


Pto + S- & o e^o-t) 


t > to 


(18) 


a — 5 T 


(3 t0 5 e 5 6o-0 


t>t Q 


(19) 


(3 t0 + S- (3 t0 e^o-t) 

When <5 = 0 the number of moles of carbon and metal oxide are in stoichiometric balance. The 
above results may be specialized for <5 = 0 to 


(3 = a = 


Pt 0 


1 + Pto(t — to) 


t>t 0 


( 20 ) 


It can be shown that as f ^ oo 


/5 — > 0, a — > 5 for <5 > 0 


( 21 ) 


P — > —5, a — > 0 for <5 < 0 


( 22 ) 


P — > 0, a — > 0 for 5 = 0 


(23) 
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FIGURE 4: Plot of the scaled concentrations of carbon (a) and metal oxide ((3) in the melt and 
the scaled rate of production of carbon monoxide (a/3) for various values of A, with carbon loss 
control (solid curves), for a case with <5 = 0. Methane flow is turned off beyond the time at which 
the curves for a and (3 intersect. Dashed curves correspond to a case where the methane flow is 
not turned off (same as in Figure 2). 
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(a) (b) 

FIGURE 5: (a) the scaled time t c at which the carbon monoxide production rate is a maximum and 
(b) the corresponding scaled maximum rate of production of CO (a(t c )(3(t c )), both as a function 
of A, with carbon loss control employed. The dashed curve corresponds to a case where methane 
flow is not turned off (same as in Figure 3). 


Therefore, when 5 < 0, there is potential for zero carbon loss in the melt, if the chemical conver- 
sion is carried out for a sufficiently long time. 


Plots of the scaled concentrations of carbon, metal oxide, and the rate of production of carbon 
monoxide are shown by the solid curves in Figure 4, where carbon loss control is implemented 
with <5 = 0. Note that <5 = 0 a to = (3 to . Thus, the methane flow is turned off beyond the time 
at which the curves for a and (3 intersect. For the dashed curves in Figure 4, the methane flow is 
not turned off. These are the same as the curves shown in Figure 2. Comparing the solid and the 
dashed curves, we see that with carbon loss control, the amount of carbon remaining in the melt (a) 
is drastically reduced after the methane flow is turned off. As a consequence, the corresponding 
metal oxide concentration (/3) is larger than that without control of carbon loss, and it takes a 
longer time for the conversion to be completed. When A < 1, however, the rate of production of 
carbon monoxide in general (proportional to a/3), and in particular the peak CO production rate, 
are not significantly affected. In contrast, when A = 10, the carbon monoxide production rate 
is significantly reduced when carbon loss control in implemented. The dependence of the time 
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t c at which the CO production rate is a maximum, and the scaled maximum CO production rate 
a(t c )/3(t c ) on A are shown in Figure 5 when carbon loss control is implemented. In this figure, t c 
satisfies Equation (13) for A < 0.842, and equals j for A > 0.842. 

3.3 Growing melt 

We now analyze a case where the melt volume increases with time. Specifically, we will permit 
the volume to increase as the square root of time. This would be the case when the regolith (of un- 
bounded extent), initially at a constant temperature below its melting temperature, is heated so that 
its surface temperature is instantaneously raised and maintained above the melting temperature. 
The rate of growth of the melt volume will depend on the latent heat and other thermal properties 
of the regolith. In our analysis, we will assume that the melt volume V(t*) = V 0 yj^, where V 0 
and t 0 are known reference values. We will assume that the surface area A of the melt surface is 
a constant, which is valid when the heat transfer is predominantly one-dimensional. We retain all 
other assumptions for the chemical process discussed before. 


The conservation of mass yields the following equations for the molar concentration of carbon 
and metal oxide in the melt. 

j-(VFa*) = AjC ™' /r ° - kV¥a‘ 0 * (24) 

at* Vn 


- a kVt* a* (3* (25) 

We will permit a time delay between the start of melt growth (t* = 0) and the initiation of methane 
flow (t* = ti). The initial conditions are a* = 0 and (3* = /3 0 at t* = f*. Dimensionless variables 
are defined as follows. 


a 


a = 


/ 3 * 


a , > P = P= (2 kp 0 t*) 1/2 , Pi = ( 2kf3 0 t i ) 1/ 2 

Po /a (3 0 


The governing equations in terms of these variables are 


d 

dp 


(P a) 


(2k/3 0 t 0 ) 1 ^ 2 X p — p 2 a (3 


(26) 


(27) 
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^(p/3) = 1 -p 2 a(3 


(28) 


a(pi) = 0, Pfa) = 1 (29) 

Equations (27) and (28) may be partially solved and written as 

y 

pa = pP+ — (p 2 -p^) -p (30) 

^(p/3) = 1 ~pP (pP + j(p 2 ~P 2 i) -pj ( 31 ) 

where A' = (2k/3 0 t 0 ) 1 ^ 2 X. A is the same parameter defined earlier (with V 0 used for the melt vol- 
ume). The dimensionless parameter k/3 0 t 0 is the ratio of the time scale for the growth of the melt 
volume to the time scale for the reaction in the melt. 


The instantaneous CO production rate is 


dXco 

dt* 


V k a* (3* 


Vokffi n 

a(2kf3 0 t 0 y/ 2 Pa 


(32) 


Thus, the CO production rate is proportional to p a /3. Equation (31) for p f3 is a form of the 
Riccati equation, but does not appear to have a closed-form solution. We have solved Equation 
(31) numerically. Plots of the scaled concentrations of metal oxide (/3), carbon (a), and the rate of 
production of CO (p a (5) are shown in Figure 6 for various values of A' when the methane delay 
time ti = 0. A log-log scale has been used in these plots to concisely capture the various features 
of the curves shown. Many of the qualitative features of the curves in Figure 6 are the same as 
those in Figure 2 for a constant volume melt. The metal oxide concentration decreases monotoni- 
cally with time. For large times, the carbon concentration increases with time because there is no 
carbon consumption once the metal oxide is depleted. The rate of production of CO is initially 
zero, attains a maximum at some time, and finally falls off when the metal oxide is depleted. An 
interesting feature in Figure 6 not present in Figure 2 is the behavior of carbon concentration with 
time for small values of A', i.e., A' = 0.01 and 0.1. The carbon concentration has a local maximum 
and a local minimum before it increases monotonically for large values of time. This behavior 
may be explained as follows. Initially the carbon concentration a is zero. Deposition at the surface 
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FIGURE 6: Case of a growing melt - Plot of the scaled concentrations of carbon (a) and metal 
oxide (/3) in the melt , and the scaled rate of production of carbon monoxide (p a (5) for various 
values of A', versus p (proportional to square root of time, Eq (26)). 
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tends to increase a, while consumption due to the reaction and dilution due to increasing melt 
volume both tend to decrease a. a would tend to increase linearly with time for a constant rate of 
deposition without consumption or dilution. For small time, the reaction is slow, and the constant 
rate of carbon deposition overpowers its dilution, as the melt volume increases at a slower rate 
(growing as y / f). As time increases, the carbon consumption due to the reaction increases. Thus, a 
attains a maximum and subsequently decreases for a period of time until the reaction slows down 
due to metal oxide depletion, a thus attains a minimum. Thereafter, carbon deposition once again 
overpowers dilution, and a increases with time. Between the maximum and minimum of a. Figure 
6 shows that the rate of production of CO is practically constant. 

When A' is sufficiently large (A' > 1), consumption and dilution of carbon are both weak com- 
pared to its deposition, and no local extrema in a are present in Figure 6 . Control of carbon loss 
in the melt, especially for large A', may be implemented along the lines discussed in the previous 
section, if desired. 

We explored the effect of delaying the initiation of methane flow, and the results are shown in 
Figure 7 for A' = 0.01 and 0.1. The maximum carbon concentration achieved after the methane 
flow is turned on is lowered. The longer the delay in methane flow initiation, the lower is the peak 
of carbon concentration. Analogous to the carbon loss strategy discussed earlier, manipulating the 
methane flow initiation time can be used as a method of controlling carbon build-up in the melt 
and possible carbon crust formation. The peak rate of production of CO is not affected by delaying 
the initiation of methane flow. 

In the models that we have described above for both constant and growing melt volume, the rate 
constants k and 7 are crucial physical properties that the model requires as input. These rate con- 
stants are not readily available in the literature and must be obtained by comparison of the model 
predictions with experimental results. Results from well-controlled experiments, however, do not 
seem to be available. Preliminary estimates k = 4.5 x 10~ 8 ^ and 7 = 0.04 — have been ob- 
tained for a lunar regolith simulant (see Balasubramaniam et al., 2009 for details) that best fit the 
test data obtained by ORBITEC. In these tests, however, the flow rate of methane in the gas phase 
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(a) A' = 0.01 




(b) A' = 0.1 

FIGURE 7: Effect of delay of methane flow initiation for a growing melt - Plot of the scaled 
concentrations of carbon (a) and the scaled rate of production of carbon monoxide (p a (5) versus 
p for pi = 0, 0.5, 1, 1.5, 2, 2.5, for (a) A' = 0.01 and (b) A' = 0.1. For each curve, the value of Pi 
corresponds to the value of p at which the curve crosses the p axis. 
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was sufficiently low that a substantial amount of the incoming methane was pyrolyzed. Thus, the 
model assumption that gas phase concentration of methane at the surface can be assumed to be 
constant and approximated by its value at the flow inlet is questionable in these tests. These values 
for the rate constants must be used with caution and as preliminary estimates only. 

4 Conclusions 

We have developed a model which describes the chemical conversion of lunar regolith in a car- 
bothermal reactor to produce carbon monoxide, and eventually oxygen. Methane flows over a 
molten regolith zone in the reactor continuously at a fixed temperature and pressure to maintain its 
concentration. Our model assumes a two-step process for the chemical conversion. First, methane 
is pyrolyzed as it comes in contact with the hot, molten regolith surface. The solid carbon formed 
after pyrolysis is deposited on the melt, and mixes with it uniformly. The second step is the reac- 
tion which takes place between carbon and the metal oxide in the melt to yield carbon monoxide 
that bubbles out of the melt. We considered both constant volume and growing melt cases. For 
a constant volume melt, our model shows that a single parameter captures the effects of melt size 
relative to its surface area, reaction rate constant, pyrolysis rate constant, and the concentrations of 
methane and metal oxide on the rate at which chemical conversion takes place in the reactor. For 
a growing melt, the ratio of the time scale for melt growth to the chemical time scale also appears 
in the dimensionless parameter. Simple strategies to control the amount of unreacted carbon in the 
melt by turning the flow of methane off at a suitably chosen time, and to control the build up of 
carbon at the beginning of the process when the melt is growing are discussed. Rigorous validation 
of the model predictions with experimental results needs to be performed. 
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